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SmiLAIION  M3DEL 


A. .  INTRODUCTION 

This  appendix  documents  the  mathematical  model  used  in  the  unpiloted 
(digital)  and  piloted  high  AOA  simulations  of  a  generic  high-speed^  tvrin-Jet, 
fighter-type  airplane.  The  model  represents  a  relatively  "clean"  aircraft 
(no  flap,  slat,  or  gear  deflection)  based  upon  an  r-4j,  with  external  stores 
consisting  of  one  centerline  fuel  tank  and  one  air-to-air  missile .  Four 
different  generic  configurations  are  generated  by  changing  a  limited  number 
of  aerodynamic  coefficients. 

The  appendix  is  organized  as  follows:  Section  B  presents  the  equations 
of  motion  including  a  discussion  of  assumptions,  axis  systems  and  transfor¬ 
mations.  The  equations  for  the  aerodynamic  force  and  moment  coefficients 
are  presented  in  Section  C,  The  airplane  physical  characteristics  and  con¬ 
trol  system  are  described  in  Section  D<  Section  £1  summarized  the  data  sources. 
The  aerodynamic  coefficient  data  are  presented  in  Section  F.  The  coefficient 
"lookup"  tables  are  listed  subsequentlally  in  Table  9  and  plotted  in  Pigs.  5 
to  28.  Validation  of  the  model  with  flight  test  data  is  given  in  Section  0, 

B.  EQUATIONS  OF  MDTION 

The  nonlinear  slx-degree-of-freedom  aircraft  equations  of  motion  used  in 
the  simulation  programs  are  consistent  with  the  following  assumptions, 

1  •  ABBumptions 

a.  The  airframe  is  assumed  to  be  a  rigid  body. 

b.  The  earth  is  assumed  to  be  fixed  in  inertial  space, 

c.  The  mass  and  mass  distribution  of  the  vehicle  are 
assumed  to  be  constant. 

d.  The  aircraft  has  a  plane  of  symmetry. 
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e.  Effects  associated  with  rotation  of  the  vertical 
relative  to  inertial  space  are  assumed  negligible j 
the  magnitude  of  the  gravity  vector  is  assumed 
constant. 

2.  Axis  Systems 

All  the  axis  systems  (see  Fig.  i )  have  their  origin  at  the  aircraft 

center  of  gravity  (o.g.)«  The  Earth  Axis  system  is  oriented  with  the  Zjj 
axis  along  the  local  gravity  vector  and  the  Yg  axes  in  a  horizontal 
plane  with  arbitrary,  but  fixed,  direction. 

The  Body  Axes  ere  fixed  In  the  aircraft  with  the  axis  positive 
forward  along  the  fuselage  reference  line.  The  Zg  axis  is  in  the  plane 
of  symmetry,  positive  down,  and  the  Yg  axis  is  perpendicular  to  the  plane 
of  symmetry,  positive  out  the  right  wing,  llie  Body  Axes  are  located 
relative  to  the  Earth  Axes  by  conventional  Euler  Angles  Y,  ®,  and  <I>. 
y  is  a  rotation  of  the  Body  Axes  from  the  Barth  Axes  about  the  Zg  axis; 
e  Is  a  rotation  about  the  Intermediate  Yg  axis;  and,  finally,  is  about 
the  axis.  All  angles  and  anguD.ar  rates  are  positive  In  a  righthanded 
sense. 

The  final  axis  system  is  called  Flight  Path  Axes.  They  are  aligned 
with  the  aircraft  Inertial  velocity.  For  the  still-air  case  considered 
here  they  are  identical  to  Wind  Axes  (usually  aligned  with  the  aircraft 
velocity  relative  to  the  air  mass).  Herein  we  use  the  more  conventional 
W  subscript  and  Flight  Path,  Flight,  and  Wind  Interchangeably.  The  axis 
lias  along  the  aircraft  total  velocity  vector,  Vj;  the  Z^  axis  is  in  the 
aircraft  piano  of  symmetry;  and  the  Y^  axis  completes  the  right-handed 
orthogonal  set.  The  Wind  or  Flight  Axes  are  located  relative  to  the  body 
by  the  angles  of  attack,  a,  and  sideslip,  B;  a  Is  a  rotation  about  the 
axis  and  p  a  rotation  about  tV.e  Zyf  axis  (if  a  =  o) . 

The  transformation  from  Flight  Path  Axes  to  Body  Axes  is  given  by; 


cos  a  cos  B 

-cos  a  sin  B 

-sin  a 

% 

a 

sin  B 

cos  B 

0 

2b 

sin  a  cos  B 

-sin  a  sin  B 

cos  a 

• 

. 

2 


■x 

. . . . .  ’  ■ 


3 .  Moment  Equations 


P,  and  R  are  the  instantaneous  components  of  the  aircraft  total 
inertial  angular  velocity  vector  in  the  Xg,  Yg,  Zq  Body  Axes.  The  deriva¬ 
tives  of  these  components  with  respect  to  time  are  related  to  the  applied 
moments  about  the  center  of  mass  by: 


where 
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Ijj,  ly,  Igi  and  a*"®  "the  moments  and  product  of  inertia  with  respect  to 
the  Body  Axes.  X ,  ifl,  and  7L  are  moments  about  the  Body  Axes  due  to 
aerodynamics  and  aircraft  thrust. 


If.  Force  Equations 

Flight  Path  Axes  force  equations  are  used  for  computational  efficiency. 
The  derivatives  of  a,  p,  and  the  aircraft  total  inertial  translational  veloc 
ity,  V^,  are  related  to  the  external  forces  by: 


a  =  Q  ~  tan  p(P  cos  a  +  P,  sin  a)  +  Z^/(mVip  cos  P) 

P  =  P  sin  a  -  R  cos  a  +  Y^/(aV,j,)  (Vi 

Vj  "•  3^/® 

and  are  the  components  of  the  total  external  forces  (aerodynamic, 
thrust,  and  gravity)  along  the  Flight  Path  Axes*  in  is  the  aircraft  mass* 

The  a  equation  is  limited  in  accuracy  to  small  sideslip  angles.  The 
large  excursions  in  angles  of  attack,  and  sideslip  encountered  in  the  piloted 
simulation  necessitated  a  change  in  definition  of  a.  The  Body-Axis  equation 
for  vertical  acceleration,  w,  was  substituted: 

W'  *=  QV^  cos  p  -  Vqi  sin  p  (P  cos  a  +  R  sin  a)  +  2i^/m. 

and  the  definitions  for  a  and  p  were  a  «  sin"^  'w/V^,  p  =  sin“^  v/Vj* 
These  excursions  were  not  experienced  in  the  preliminary  analytical  work  per^ 
formed  in  preparation  for  the  piloted  simulation. 

S .  Thrust  Geometry 


The  thrust  is  assumed  to  lie  in  the  aircraft  plane  of  symmetry  and  is 
oriented  with  respect  to  the  Body  Axes  by  the  thruet  inclination,  and 
offset,  Zy  as  indicated  in  the  sketch  below. 


X| 


6»  Body  Axis  Moments 

The  roll, ,2^  ,  and  yaw,7L,  moments  are  due  to  aerodynamics  only,  while 
the  pitching  moment, 77Z,  Includes  a  thrust  term.  These  are  given  by; 


=  qSbC^  ;  ^ 


qScCm  +  ZjT 


)  ^  a  qSbCf, 


(5) 


- :::: . . . . 


Tl  is  the  dynamic  pressure;  b,  5,  and  S  are  the  reference  span,  chord,  and 
wing  area,  respectively.  C^,  C^,  and  Cj^  are  the  total  body-axis  non- 
dimensional  aerodynamic  moment  coefficients  (referenced  to  the  c.m.). 

7>  Flight  Path  Axis  Forces 


The  forces  along  the  Plight  Path  Axes  have  components  due  to  thrust  (T), 

weight  (mg),  and  aerodynamics.  They  are  given  by; 

,  • 


% 


T  cos  p  cos  (a  +  4q) 

+  mg|aoa  0  cos  4  sin  a  cos  P  -  sin  0  cos  a  cos  P  +  cos  0  sin  $  sin  6 
+  sin  p  -  D  cos  p 


-T  sin  P  cos  ( a  +  Jg) 

+  mgjooa  ©  sin  ®  cos  P  +  sin  0  cos  a  sin  p 
+  cos  p  +  D  sin  P 


-  cos  0  cos  sin  a  sin  p 


Zyj  ■  -T  Sin  (a  +  5g)  +  mg^sln  0  sin  a  +  cos  e  cos  ®  cos  a>  -  l 


L,  D,  and  Y^  are  the  aerodynamic  forces  and  are  given  by; 

L  =  ^SCl 

D  =  -gscD  (7) 

Ya  -  gSCy 

Cj^,  Cjj,  and  Cy  are  the  lift,  drag,  and  side  force  coefficients.  Lift,  L, 

and  drag,  D,  are  assumed  in  the  plane  of  symmetry;  lift  positive  up  and  per¬ 
pendicular  to  the  velocity  vector,  drag  positive  aft  and  along  the  projec¬ 
tion  of  the  total  velocity  In  the  plane  of  symmetry.  The  side  force  Y/^  iti 
that  component  of  the  total  aerodynamic  force  which  is  perpendicular  to  the 
plane  of  symmetry;  it  lies  along  the  Yg  axis  and  is  positive  out  the  right 
wing. 
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It  should  be  noted  that  the  drag  and  side  force  are  not  the  same  as  the 
X  and  Y  Flight  Path  Axes  forces  for  non-zero  sideslip^  see  Fig.  2.  The  aero¬ 
dynamic  forces  instead  have  been  defined  in  terms  of  the  axes  in  which  aero¬ 
dynamic  data  are  moat  frequently  available. 

8t  Body  Axis  Euler  Angle  Rates 

The  Body  Axis  Euler  angle  rates  are  related  to  the  angular  velocity 
components  by; 


f  ■  (R  eos  0  +  Ci  sin  0)/cos  © 

0  ■  Q,  008  0  -  R  sin  0  (8) 

♦  ■  P  +  T  aln  0 

9.  Body  Axis  Velocities 

The  componemts  of  the  total  velocity  along  the  Body  Axes  are: 

U  -  cos  a  oos  B 

V  -  Vj  sin  B  (9) 

W  ■  Vj  sin  a  cos  B 

10.  Earth  Axis  Velocities 

The  Body  Axis  velocities  are  transformed  through  the  Euler  angles  to 
yield  the  Earth  Axis  velocities. 

« 

Xg  ■  U  cos  0  cos  ▼  +  W(cos  0 

+  V(Bin  0 

Yjj,  ■  U  cos  0  sin  T  +  V(8in  0 

+  W(co8  ♦ 

Zjjj  ■  -U  sin  0  +  V  sin  0  cos  0  +  V  cos  <S  eos  0 
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sin  0  cos  V  +  sin  0  sin  v) 

sin  0  cos  T  —  cos  0  sin  v) 

sin  0  sin  f  +  cos  0  cos  v)  (IO) 

sin  0  sin  T  —  sin  0  cos  v) 


Figure  g.  Relation  of  Aerodynamic  Force  Axis  Syat 

to  Wind  Axes 


Ill  Aircraft  Controls 

The  aircraft  is  assumed  to  have  pitch,  roll,  and  yaw  controls,  i.e,, 
®stab'  ^spi  ^r*  Positive  deflections  and  travel  limits  for  these 
controls  are  given  in  Section  D  of  this  appendix. 

12.  Auxiliary  Variables 

Additional  variables  which  are  useful  in  aircraft  dynamics  studies 
are  given  below. 

a.  Flight  Path  Inclination  Angle,  7 


tan  r  -  +  y2 

sin  7  “ 


sin  7  sin  8  cos  a  cos  p 
•-cos  8  sin  4  sin  ^ 

-cos  8  cos  4  sin  a  cos  p 


b.  Altitude  Rate,  K 


H  -  -Z 


c.  Body-Mounted  Accelerometer  Signals,  at  Center  of  Gravity 

axcgj  ®'ycg>  ^Zcg 

r  av  n  r®®”  ®®®  ^  —cos  a  sin  p  -sin  a  "1 


sin  p 


cos  p 


sin  a  cos  p  -sin  a  sin  p  cos  a 


X  Vt(P  -  P  sin  a  +  R  cos  a) 

Vj  cos  p[A  -  Q  f  tan  P(P  cos  a  +  R  sin  a)] 


-sin  8 


-  g|  cos  6  sin  4 


COB  6  cos  4 


or 


T  cos  Iq  —  D  cos  a  +  L  sin  a 

Ya 

-T  sin  Sq  -  D  sin  a  -  L  cos  a 


d.  Body- Mounted  Accelerometer  Signals  at  Arbitrary  Point 


ijj,  ly,  and  ^2  coordinates  of  the  arbitrary  point  in  Body  Axes, 

e.  Stability -Axis  Angular  Velocities, 

A  body-fixed  Stability  Axes  system  (Ref.  l)  may  be  defined  which  is 
located  with  respect  to  Body  Axes  by  the  trim  angle  of  attack,  a^.  The 
components  of  the  aircraft's  total  inertial  angular  velocity  in  body-fixed 
Stability  Axes  are  given  by: 

Pg  a  P  cos  ttQ  +  R  sin  ttg 

Qs  “  (15) 

R3  ■  -P  sin  ag  +  R  cos 


f.  Stability-Axis  Euler  Angles,  Yg,  ©g,  <I>g 

The  body-fixed  Stability  Axes  system  defined  above  can  be  located  relative 
to  Earth  Axes  by  a  conventional  Euler  Angle  set,  Yg,  ©g,  The  Stability 

Axis  Euler  Angles  are  related  to  the  Body  Axes  Euler  Angles  and  the  trim 
angle  of  attack  by: 


tan  Y. 


cos  ttp  sin  Y  cos  0  +  sin  ao(sin  Y  sin  0  cos  4  —  cos  Y  sin  t) 
cos  ttg  cos  Y  cos  0  +  sin  a^cos  Y  sin  ©  cos  *  +  sin  Y  sin  4) 


sin  0u  =  cos  tto  sin  ©  -  sin  a-  cos  ©  cos  0 


tan  0g  -  oos 


cos  Oq  sin  ©  —  sin  cos  ©  cos 
cos  aQ  cos  0  cos  <!>  +  sin  sin  © 


cos  ©  sin 


tan 


3  cos  ttg  cos  ©cos  ®  +  sin  Oq  sin 


g.  Stability-Axis  Euler  Angle  Rates,  Ys,  0g, 


tg  “  (Rg  cos  ®g  +  Qg  sin  ®s)/cos  Sg 
(Bg  ■  Qg  cos  <I>a  -  Rg  sin  «a 
®8  ®  ^s  '^’^s  sin  0g 

h.  Flight  Path  Azimuth  Angle,  X 
tan  X  =  Y/X 


C.  AERODYNAMIC  FORCE  AND  MOMENT  COEFFICIENT  EQUATIONS 


Aerodynamic  forces  and  moments  have  been  defined  in  terms  of  the  total 
non-dlmenaional  aerodyn.'jraic  coefficients,  i.c.,  the  stability-axit,  force 
coefficienta  Cl,  C^,  Cy  and  the  body-axis  moment  coefficients  C^i. 

As  sho-wn  in  Fig.  3,  these  are  functions  of  the  aircraft  state,  the  coni rol 
surfaces,  tabular  coefficient  and  derivative  data,  and  con;-. bant  parameters. 
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Coefficients  and  Derivatives 


The  coefficient  equations  are  given  in  Table  1 .  A  summary  of  the  tabu¬ 
lar  functions  is  given  in  Table  2  (the  data  is  given  in  Section  F)  and  the 
constant  parameters  are  presented  in  Tables  5  and  5.  As  indicated  in  Table  2, 
all  tabular  data  are  given  for  ^-deg  increments  in  a  for  either  0  deg  <  a  < 
no  deg  or  0  deg  <  a  ^  45  deg.  Linear  interpolation  between  adjacent  data 
points  was  used  to  obtain  derivative  values  for  intermediate  values  of  a. 

All  coefficients  and  derivatives  except  ClbasiC  ^®BASIC  symmetric 
functions  of  a  and  the  absolute  value  of  a  was  used  to  obtain  the  derivative 
value.  is  an  anti- symmetric  function  of  a;  as  indicated  in  Table 

|a|  is  used  to  obtain  the  table  value  and  the  sign  of  the  appropriate  term 
is  changed  for  negative  a'  a .  In  general  for  |  a] ' s  which  exceed  the  range 
of  the  tables,  the  derivative  value  for  the  maximum  ia|  given  in  the  table 
is  used.  This  is  shown  explicitly  in  Table  1,  where  the  complete  equations 
have  been  given  for  each  a  range.  It  should  be  noted  that,  except  for  Cp,, 
the  equations  themselves  are  not  dependent  on  the  a  range.  is  not 

quite  anti- symmetric}  the  explicit  equations  for  lift  coefficient  at  nega¬ 
tive  angles  of  attack  are  given  in  Table  1 , 

is  the  only  two-dimensional  function  used;  the  table  contains  values 
for  I p1  from  0  deg  to  50  deg  in  5  ^eg  increments.  For  Ipj  >  30  deg,  the  table 
values  for  p  >=  io  deg  were  used.  The  organization  of  the  data  is  discussed 
further  in  Section  F. 

As  Table  2  indicates,  there  are  two  sets  of  data  for  each  of  bhe  deri¬ 
vatives  C_0p^,  C^pj,  and  Selecting  the  proper  combinat5.ons  of 

these  derivatives  results  in  the  desired  four  basic  aircraft  configurations 
as  indicated  in  Table  4, 
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T-AHTE  1  .  AEROBYKAMIC  COEFFICIEHT  EQUAIIONS 
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TABLE  1 .  AERODTOAMIC  COEPB’ICIEIW?  EQUATIONS  (Concluded) 


b)  Body  Axis  Moment  Coefficients 
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TABLE  2.  0RGAKI2ATI0N  OF 
TABULAE  AERODYNAMIC  DATA 


PARAHETBR3 

FUNCTION 

OP 

uNm 

MNEMONIC 

AND  ARRAY 
DIMENSIONS 

LINE 

NUMBERS 

a 

QfllD** 

°i.basic 

a 

- 

CIiBAS(23) 

2-8 

o.(5.)llo. 

°^6aTAB 

a 

l/dag 

0LSTAB(85) 

9-15 

o.(5.)llo. 

°J>BASIC 

a 

- 

0DBAS(S3) 

16-22 

0.(5. )no. 

% 

a 

l/d«g 

CYB(25) 

25-29 

0.(5. )110. 

®yBr 

a 

l/d«g 

CYDR(lO) 

30-33 

0.(5. )45. 

a 

l/d«g 

CRB1(S3) 

34.40 

0.(5. )110. 

Olpg 

a 

l/d«g 

CRBS(23) 

41-47 

0.(5.)110. 

°^P1 

a 

l/rad 

CRPI(IO) 

48-51 

o.(5.)45. 

°ip2 

a 

l/r»d 

CRI>a(lO) 

52-55 

0.(5. )45. 

Qb 

l/rod 

CRR(lO) 

56-59 

o.(5.)45. 

Q 

l/dag 

OHDA(lO) 

60-63 

0.(5.)45. 

a 

l/dag 

CRDSP(lO) 

64-67 

o.(5.)45. 

Ciftr 

a 

1/dag 

crdr(io) 

68-71 

0.(5.145. 

°«BASIC 

a 

- 

CMBAS(23) 

72-78 

0.(5. )110. 

ACni 

a,  e 

- 

DCM1(10,7) 

79-95 

0.(5. )45. 

AOng 

a,  p 

M 

DCM2(10,7) 

96-112 

0.(5.145. 

Qniq 

a 

l/rad 

OM^ClO) 

115-116 

0.(5.145. 

Offli 

a 

l/r«d 

cmad(io) 

117-120 

0.(5. 145. 

a 

l/deg 

CMSTAB(25) 

121-127 

o.(5.)llo. 

a 

l/dag 

cmda(io) 

126-151 

0.(5.145. 

a 

l/deg 

CMD3P(10) 

152-135 

o.(5.)45. 

Cnp, 

a 

1/dag 

CNB1{10) 

136-139 

0.(5.145. 

°"P8 

a 

l/deg 

CNB2(l0) 

i4o-i43 

0.(5.145. 

®np 

a 

1/rad 

CHP(IO) 

1 44-1 47 

0.(5.)45. 

On, 

a 

l/rad 

CNH(lO) 

148-151 

0.(5.145. 

a 

l/deg 

CNDA(  1 0) 

152-155 

0.(5.145. 

Oflp 

a 

l/dag 

CNDSP(IO) 

156-159 

0.(5.145. 

°r 

a 

l/deg 

cndr(io) 

1 60- 1 65 

0.(5.145. 

*  0.(1;. )nOi  irdlcutud  valuau  nr«  Klvan  for  a  from  0  don  t.u  110  dop, 

In  5  dug  int;roim!rit,B,  l.o.,  0",,  irP.,  .  100".,  in*.”.,  llo". 
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TABLE  3.  CONSTANT  AERODYNAMIC  COEFFICIENTS 


^'^DgTORES^“^ 

.0037 

Cyg^(l/deg) 

-.0001 67 

Cy_  (l/deg) 

^sp 

-.00006 

K°“^Bstab(i/eieg2) 

.00011 

TABLE  4.  AIRCRAFT  CONFIGURATTaN/DERmTIVE  COMBINATIONS 


AIRCRAFT 

DERIVATIVE  INDEX 

CONFIGURATION 

i 

i 

ACmi^ 

k 

1 

A 

1 

1 

1 

1 

B 

1 

2 

1 

1 

C 

2 

1 

1 

1 

D 

1 

1 

2 

2 

S.  AStlLAM  PHYSICAL  CHARACIEBISTICB  ANB 
COmOL  SYSTEM 

1 .  Physical  Chaxacteristlcs 

This  section  describes  the  pertinent  dimensions  of  the  airplane  (weighty 
wing  area,  inertias,  etc.)  and  defines  the  control  surfaces.  The  model  is 
based  upon  the  F-4-J,  but  is  considered  to  be  representative  of  a  generic 
fighter-type  airplane . 

The  major  dimensional  parameters  for  this  model  are  listed  in  Table  5, 
the  control  surfaces  are  outlined  below  and  control  surface  limits  are  given 
in  Table  6. 

2.  Control  System 

Longitudinal  control  is  provided  by  an  all-moving  horizontal  tail.  Lat¬ 
eral  control  is  provided  by  a  combination  of  spoilers  and  ailerons.  The  ail¬ 
erons  deflect  downward  only;  the  spoilers  deflect  upward  only.  The  left  ail¬ 
eron  and  right  spoiler  operate  simultaneously,  as  do  the  right  aileron  and 
left  spoiler.  Spoiler  and  aileron  deflection  are  combined  by  a  simple  relation; 

6sp  “  1.^55&a 

Directional  control  Is  provided  by  a  conventional  rudder. 

The  longitudinal  control  system  is  outlined  in  Table  T,  and  the  lateral/ 
directional  control  system  is  shown  in  Fig. 


TABLE  5.  AIRPLAKE  DIMENSIONAL  PARAMETERS 


SYMBOL 

DEFINITION 

VALUE 

Ts 

Wing  span,  ft 

58.67 

0 

Wing  mean  aerodynamic  chord,  ft 

16. 04 

S 

Wing  area,  ft^ 

530 

W 

Airplane  weight,  lb 

37,000 

^ref 

Reference  center  of  gravity, 

31 

^og 

Actual  center  of  gravity,  ^0 

29.3 

Ix 

Moment  of  inertia  about  X-body  axis, 

BlUg-ft^ 

25,850. 

ly 

Moment  of  inertia  about  Y-body  axis, 

slug-ft^ 

127,1^00. 

Moment  of  inertia  about  Z-body  axis, 

alug-ft2 

iij-6,000. 

Ixz 

Product  of  inertia  about  X-Z  axes, 

slug-ft2 

2210. 

Thrust  inclination,  dog 

5*25 

Thrust  offset,  ft 

-0.336 

H 

Reference  altitude,  ft 

15,000 

19 


TABLE  6. 

CONT'ROL  SURFACE  DEFINITIONS 


SURFACE  DEFIiECTION  LIMITS,  DEG 


^staT3>  Positive  TED' 
Positive  Left  TED 
Bjjp,  Positive  Right  TEU 
Positive  TEL 


21  up,  9  <iown 
0  up,  30  dovn 
43  up,  0  down 
130 


TABLE  7. 

LONGITUDINAL  CONTROL  SYSTEM 


•  No  SAS  or  CAS 

•  Basic  F-4  manual  control  system  gain, 

*^stab  ^long  =  '+•■'6  deg/in, 

•  Feel  system  —  simple  spring/damper 
(gradient  *''3.5  Ib/in) 


Breakout  —  1  .5  to  2  lb 
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S.  BAIA  SOUROSS 


This  section  outlines  the  sources  for  the  aerodynamic  data  presented  in 
this  appendix.  Reference  l6  was  used  as  the  initial  data  base  for  this  simu¬ 
lation  model.  The  Ref,  l6  data  covers  an  a  range  of  -5  deg  to  3l  deg  and  a  p 
range  of  -10  deg  to  10  deg.  Wind  tunnel  tests  of  an  F-h  model  (Ref,  2),  cov¬ 
ering  an  a  range  of  ~10  deg  to  110  deg  and  a  p  range  of  -4o  deg  to  deg^ 
were  the  primary  sources  for  the  extension  of  static  coefficients.  Dynamic 
coefficients  were  extrapolated  based  on  trends  of  data  available  Itora  other 
sources  to  60  deg  angle-of-attack,  with  the  intent  of  making  such  coefficients 
simple  to  mechaniEe. 

Because  of  the  nature  of  the  reference  (an  STI  working  paper),  some  dis¬ 
cussion  is  in  order  on  the  data  which  was  used  to  develop  the  initial  data 
set.  Three  separate  sources  were  utilized  in  Ref,  16. 

The  first  two  seta  of  data  were  supplied  to  STI  for  a  previous  study. 

One,  Ref,  3,  was  based  on  the  spin  study  data  (Ref.  1+) .  The  second  was  sup¬ 
plied  by  NASA  Langley  and  was  based  on  several  NASA  tunnel  investigations 
(e.g,,  Refs.  5  and  6).  The  Ref,  3  tabulation  was  a  function  of  a,  p,  and 
Sgtab  range: 

0  <  a  <  90  deg  )  0  <  p  <  4-0  deg  ;  -21  <  63^.3, <  0  deg 
The  NASA  tabulation  was  a  function  of  a  and  p  for  the  range: 

-10  <  a  <  120  deg  j  -H-0  <  P  <  +4o  deg 

Both  sets  are  limited  to  the  very  low  speed  regime  (M  »  0.2).  As  a  result 
of  comparison,  several  erroneous  data  points  were  discovered  in  the  Ref.  3 
data  and  were  reported  to  the  AFFDL. 

The  third  set  of  data  (Ref.  7)  was  obtained  from  the  NASA  Langley  Differ¬ 
ential  Maneuvering  Simulation  (DMS)  Investigation  of  the  F-4j  and  slatted  F-il-E 
aircraft.  Several  key  damping  derivatives  were  updated  from  Ref,  7  based  upon 
NASA  correlation  between  several  sets  of  tunnel  data  and  evaluation  of  aircraft 
response  obtained  in  the  simulation.  These  data  also  are  in  look-up  table  for¬ 
mat  with  coefficients  as  functions  of  a,  M,  and  H  over  the  range: 

-5  <  a  <  30  deg  }  .2  <  M  c  2.4  ;  15,000  <  H  <  45,000  ft 

The  sideslip  coefficients  ere  valid  over  the  range  p  =  ±4o  deg.  Thrust  and 
drag  effects  are  also  modeled  in  detail. 
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As  originally  received,  the  data  did  not  coalesce  at  0.2  M,  The  data 
of  Ref,  4  vere  therefore  given  the  itiDst  weight  at  this  low  Mach  and  were 
smoothed  into  the  DMS  data  by  about  0.4  M.  The  resulting  data  provide  an 
excellent  model  of  the  aircraft  over  the  range: 

-4  <  a  <  30  deg  j  -10  <  p  <  10  deg  j  0.2  <  M  <  1 ,0  ;  0  <  H  <  25,000  ft 

Besides  extension  of  coefficients  for  high  angles  of  attack  and  sideslip, 
several  changes  were  made  to  the  Ref.  l6  data  to  simplify  use  of  the  data. 
These  included  normalizing  the  data  for  trim  Mach  Number  (determined  from 
previous  digital  simulations  using  the  Ref.  l6  data  set),  using  a  single  alti¬ 
tude  (l5>000  ft),  and  defining  values  for  the  second  set  of  the  parameters 
Cnp,  and  ACm. 

r«  misDvxma  oomioim  data 

The  aerodynamic  coefficients  are  tabulated  in  Table  9  on  pages  25  to  27 
and  plotted  in  Pigs,  5  to  28,  All  coefficients  are  functions  either  of  a  or 
of  a  and  3,  The  aerodynamic  data  file  in  Table  9  includes  a  one-line  iden¬ 
tifier  for  each  coefficient,  plus  one  or  two  lines  describing  the  independent 
variable(s),  The  formats  of  these  identifying  lines  are  given  in  Table  8, 

All  coefficients  are  presented,  five  to  a  line,  in  3-degree  angle-of- attack 
increments . 

For  those  coefficients  which  are  functions  of  0  <  a  <  110°,  the  values 
given  correspond  to  a  as  follows: 

(O) . (20) 


(loo)  (105)  (110) 

For  0  i.  a  A  45°  the  values  are  as  follows: 


(0) 

(25) 


(20) 

(45) 


TABLE  8. 

DATA  FILE  IDENTIFICATION  LINES 


CoBfflotent  identi¬ 
fier! 

VARIABLE 

NAME 

VARIABLE 

IffllTS 

NUMBER  OF 
INDEPENDENT 
VARIABLES 

DATA  POINTS 
PER  LINE 

NUMBER 

OF  LINES 

TOTAL  NO- 
OF  POINTS 

(1  line) 

(AS) 

(AB) 

(la) 

(la) 

(la) 

(15) 

Independent  vari¬ 
able  identifier) b) 1 

INEEPKHDENI 

VARIABLE 

NAME 

IHLEPENDENT 

VARIABLE 

UNITS 

MINIMUM 

VAIAJE 

INCREMENT 

MAXIMUM 

VALUE 

NO.  OF 
VALUES 

(l-2  lines) 

(A8) 

(AS) 

(B14.6) 

(Ell*. 6) 

(EiU.6) 

(12) 

The  coefficients  DCMl  and  DCM2  are  functions  of  both  a  and  p.  In 
this  case,  a  is  incremented  first,  then  p: 

(0,0) . (20,0) 

(25,0) . (1)5,0) 

(0,5) . (20,5) 


(25,30) . (a,p) 

As  an  exeutple  in  interpreting  the  data  file,  the  file  for  Cyp  is 
given  below ! 


CYB 

PER  DEG 

1  5  5  23 

ALPHA 

DEGREES 

0., 5. ,110. ,23 

-.011 

-.011 

-.011  -.011 

-.011 

-.011 

-.011 

-.0092  -.0079 

-.0065 

-.0050 

-.0036 

-.0023  -.0008 

.0006 

.0020 

.0030 

.0030  .0030 

.0030 

.0030 

.0030 

.0030 

The  first  line  identifies  the  coafficlentj  gives  the  units;  indicates 
that  it  is  a  function  of  one  variable,  and  that  it  is  listed  five  points 
per  line  in  five  lines,  with  a  total  of  25  data  points  given-  The  second 
line  identifies  a  as  the  independent  variable;  gives  its  units;  gives  the 
minimum  value  as  0  deg,  the  Aa  increment  as  5  deg,  and  the  maximum  value 
as  110  deg,  for  a  total  of  25  separate  Aa  points- 
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TABLE  9.  AERODYNAMIC  CCEFFICIENTS 


S.O.  lOSIl  t:;<Tt:MLi|£ri  M)^  HATA.  CR£.1  it  CRIi2  CMANCit^D  IJ./20/7B 

CLBfiH  1  5  3  23 

AL.PHi'i  DEtiREKS  0  .  .  3 .  .  1.10 .  »23 

,1220  .41  ..4034  ,V053  ,9Z9i 

1.010  1.047  1.04‘?  1.024  .94'? 

.'?27  .872  .798  .477  ,344 

.41B  .:;av  .is?  .023  -.114 

-.237  -.402  -.347 


CLliiTAEt 

PGR  tiRG 

1  3  5  23 

ALPHA 

tEGREGS 

0. .5. mo. *23 

.00499 

.00494 

.00474  .00437 

.00571 

.00513 

. 00303 

.00404  ,00443 

.00438 

.00414 

.00393 

00371  .00348 

.00324 

.00300 

.00300 

.00300  .00.300 

.00300 

,00300 

.00300 

.00300 

CDBAS 

1  S  5  23 

ALPHA 

DEOREES 

0. *3. *110. .23 

.035  .044  .130  .247  .382 
.010  .455  .784  .911  J.Oll 
1.123  1.243  1.352  1.403  1.451 

1.493  1.323  1.542  1.370  1,505  - 

1.400  1.383  1.547 

CYD  PER  DEG  1  0  3  23 

AI.PHA  riEGREES  0  .  .5 .  .  110 .  .22 

“.on  -.011  -.011  -,o:u  -.on 

“.oil  -.011  -.0092  -.0079  “.0043 

“.0050  -.0034  -.0023  -.0008  .0004 

.0020  .0030  .0030  .0030  .0030 

,0030  .0030  .0030 

CYDR  F'ER  riEG  1  3  2  10 

ALPHA2  IiEQREES  0.,  5,  .45..  10 

,00151  .00142  .00134  .00130  ,00124 

,00090  .00073  .00048  .00024  0.0 

CRBl  PER  liiEfl  1  5  5  23 

ALPHA  riEGREES  0 . » 3 . 1 110 . » 23 

-.00132  -.00175  -.00234  “.00219  -.001,18 

-.00004  .00023  -.00010  -.00108  -.0010 

-.00222  -.00241  -.00255  -.00248  -.00279 

“.00291  -.00300  -.00300  -.00300  -.00300 

-.00300  -.00300  -.00300 

(:RP2  pgr  tiEC)  1  5  5  23 

ALPHA  tiGGRGES  0 .  ,S .  r  1 10 . » 23 

-,00132  -.00175  -.00234  -.00219  -.00219 

-.00219  -.00219  -.00185  -.00152  -.00180 

-.00222  -.00241  -.00253  -.00248  -.00279 

-.00291  -.00300  -.00300  -.00300  -,00300 

-.00300  -.00300  --.00300 


CRPl 

PER  RACI 

•t.  5  2  10 

ALPHA2 

OEOREGS 

0, *3. *45. .10 

-.30  4 

-.300  -. 

255  -.102 

,  100 

-.280 

-.304  -. 

250  -.164  -. 

.100 

CRP2 

PER  RAti 

15  2  10 

ALPHA2 

DEGREES 

0, *5. *45. *10 

-.304 

-.300  -. 

240  -.240 

,240 

-.240 

-.240  -. 

250  -.144  -. 

.  lOO 

CRR 

PER  RACI 

1  5  2  iO 

TABLE  9 


(CONTINUED) 


rtl.F'HALi  DHljREIi'S  0 ,  f  5 .  f  4S .  » J.(> 
.0<1,4  ,0?it  .iO£J  .170  .a:i0 
.-’(37  .3<11.  .201  ,00.'.  .OSO 


CH'OA 

PER  tiEEQ  1 

3  2  10 

Al,PHA2 

Di:nREt’8  0. 

r5, »4S. olO 

.000448 

,000439 

,000431  ,000403 

,000379 

.000334 

.000312 

.000210  .000103 

.000000 

CIRDSF' 

Pfl'R  BED  1 

3  2  10 

ALPHA2  , 

llEClIVEEa  0. 

»3.  r43,  »1,0 

.000130 

.000140 

.000093  .OOOOAA 

.  000034 

.000040 

.000030 

.000020  .000010 

.000000 

CRIIR 

PER  DliO  1 

3  2  to 

ALPHAU 

BEGREES  0, 

,<3.r43.fl0 

.000199 

.0001,43 

.000133  .00011 A 

.000092 

.000069 

.000032 

.000033  ,000017 

.000000 

CHBOS  :l  B  3  23 

ALPHA  DEOREES  0 .  f  3 .  r  1  :l 0 .  .  23 
.000  -.013  -.027  -.033  -.OBO 

-,OACI  -.003  -.144  -.126  -.242 

-.333  -.21,13  -.313  -.381  -.430 

-.47K  -.301  -.S32  -.371  -.Ali 

-,A38  -.A'W  -.729 

UCMl  2  314  70 

ALPHA2  ntiOREE-S  O..5.»43.»10 
BETA  IJIlGRtiES  0..3.i30.f  7 
.0000  ,0000  .0000  ,0000  .0000 

,0000  .0000  ,0000  .0000  .0000 


-.0080 

",0080 

■■.oouo 

-.0133 

-.0099 

.0018 

.0003 

.0000  .0000  .0000 

-.0170 

-.0170 

-.0170 

-.0233 

-.Ol'OA 

-.0204 

-.0140 

1000 

.0000  . 

0000 

-.0,312 

-.0312 

-.0312 

- . 0430 

-.0380 

-.0443 

-.0428 

-.0023 

.0000 

,0000 

-.0433 

-.0433 

-.04I5B 

-.0703 

-.0800 

-.07,48 

-.0733 

■".0087 

,0000 

.0000 

■-.0343 

- .  03,43 

■■.Or.5.43 

“.0!;i83 

-.1000 

-.0933 

-.094  4 

-.0104 

.0000 

.0000 

-.0.4  74 

-.0,474 

-.04 '74 

-.1071 

-.1200 

-.1138 

-.1133 

-.0:123 

.0000 

.  0000 

i:ic:m2  2  314  70 


ALPHA2 

BEGREES  0.»5. 

> 45, » 10 

BETA 

BEdRn- 

8  0. fS, 

>,30.  »7 

0.0  0, 

,  0  0.0 

0.0  0 

.0 

0.0  0 , 

,0  0.0 

0.0  0 

.0 

.0080 

,  . 0080 

.0080 

.0138 

.0099 

-.0018 

■  -.0005 

.  0000 

.  0000 

.0000 

,0170' 

.0170 

.0170 

.0235 

.0306 

.0204 

.01,40 

.00,10 

.0000 

.0000 

.0312 

.03:12 

.0312 

,0430 

.  0380 

.04,43 

.04251 

.oo:i3 

,0000 

,0000 

.0433 

,0433 

.  041.3 

,0703 

.0800 

.07,48 

.0733 

.0087 

.0000 

.  0000 

.  03,43 

.0343 

,  0343 

.0988 

.1000 

.0933 

.0944 

.0104 

.  0000 

.  0000 

,0,474 

.0474 

,  0474 

.1071 

.1200 

.1138 

.1133 

.0123 

.0000 

.0000 
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TABLE  9 


(CONCLUDED) 


CMC! 

PER  RAD 

1  5 

2  10 

Al.PHAi! 

degree:!! 

0 .  F  5 . 

f45. fIO 

■■•3 , 10 

-?J.41  -3 

,78 

-3.92  -3 

.loo 

”3.62 

-3.13  -3 

.13 

-3,15  -3 

.15 

CHAD 

PER  RAD 

1  5 

2  10 

ALPHA;! 

DEGREES 

0.  fS, 

f45. fIO 

“1.30 

■■■1.42  -1 

.38 

-1.63  -1 

.59 

-1 .32 

-1 , 32  -1 

.32 

“i.;3a  -1 

.33 

CMSTAB  P'l'iR  riti:o  1  ?  a  23 
<^L,F'HA  1I:iI-0RE:I“S  0,»13.»U<>.»23 
-.0:1010  -,01004  -.009131  -.00936 

-.00743  00733  •■.00677  -.00637 


-,00519  ~.0046A  -.00414  -.OOVJAl 

-.00256  -.00302  -.00150  -.00099 

0.00000  0.00000  0.00000 
CMtirt  PIER  niEW  1  5  2  10 
ALPHAa  OEORIEI-a  0.,B.»4a.il0 


-.00105  -.00092  -.00079  -.0006(3 

-.00066  -.00066  -.00066  -.00066 
CMfiSP  PER  DEG  1  B  2  10 
ALPHOa  DEGliEr.S  0 ,  ,5.  »4!!). » 10 


00828 

-.00573 

-.00309 

-.00053 


-.00066 

-.00066 


.000130  .000079  ,000035  .000007  .00000 

.000000  .000000  .000000  .000000  .000000 

CNBt  PER  OlEG  1  5  2  10 
AI.PH02  OEGRIEEG  0 .  i.  5 .  i  45 .  f  10 
.0022  ,0022  .0021  ,0013  -.0006 

-.0033  -.0040  -.0038  -.0026  -.0030 
CNlia  PER  i;iCi(3  1  5  2  10 
ALPH A2  BEGRIEES  0 .  *  M . ,  45 .  *  10 
.0032  .00:12  .0021  .0013  .0000 

-.0005  -.0009  -.0011  -.0014  -.0017 
L.'IP  PEP  P'AP  1  3  a  10 
alpha;!  riEC'IREES  0 . » 3 . » 43 .  F 1 0 
.005  -.014  .007  .oa;?  -.002 

.055  .i;22  .080  .041  .000 

CNR  PER  R'AfJ  1  5  3  10 
ALPHA2  tiEOREi:-;!:)  0,f3.»4S.»10 
•■•.37'?)  -.361  --.IJAl  -.3V0  -.502 

■-.660  -.360  -.373  -.l^i'O  .000 

CNDA  PER  DEO  1  5  ;2  10 
ALPMAl!  DEGREES  0  .  fH  .  f 4;'i .  f  10 
-.000047  -.00'00-1i;)  -.000056  -.000067  -.000068 


-.000068  -.000069  -.000070  -.000070  -.000070 

c:ni;i3F'  per  deu  i  5  a  10 

ALPHAa  DEGREE^S  0 .  f5  .  f 4'.r) .  f  10 

,00003:;!3  .0o004?33  .OO<fO3(:I0  .0000132  .0000090 

.0000044  .0000009  .0000000  .OO'OOOOO  .0000000 

(:;ni:ir  P'ER  deg  i  3  2  10 
alpha:;  degree:'.!  o.Fr..F45..io 

■-,00090  -.';>0006  -.00078  -.00073  -.00065 

-,000'.j0  -.00032  -.00033  -.00014  ,00000 


'i 
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Figure  5.  (a) 


Figure  6.  Ct  (a) 
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Figure  25.  (a) 


G.  VALmATION 


No  F-4j  flight  test  traces  of  high  AQA  are  available  to  compare  against 
dynamics  from  the  foregoing  mathematical  model.  However,  the  model  was 
originally  developed  for  use  in  a  moving-base  piloted  simulation  .  •,  'ihich 
Navy  F-4J  pilots  were  given  air  combat  maneuvering  training.  As  a  part 
of  the  checkout  and  acceptance  tests,  the  simulation  was  flown  through  var¬ 
ious  offensive  and  defensive  combat  maneuvers,  stalls,  and  departures  by 
Navy  instructor  pilots  who  indicated  it  adequately  represented  the  F-4j 
handling  and  performance.  This  provided  the  first  gross  validation  of  the 
aero  model. 

The  only  high-AOA  flight  traces  available  are  from  the  P-4lil  stall/  post- 
stall  flight  test  (Ref.  8).  McDonnell  Aircraft  Co,  indicates  that  all  hard- 
wing  models  (no  leading-edge  slats)  of  th<sj  F-4  have  approximately  the  same 
stall/departure  characteristics.  The  Air  Por'’e  (Ref.  8)  noted  the  longer 
nose  F-J+E  to  have  somewhat  less  wing  rock  tendency  than  the  C  and  D  models, 
The  F-4j  is  more  similar  to  the  C  and  D  and  therefore  migtit  also  have  more 
pronounced  wing  rook.  Also,  the  F-4E  flight  test  vehicle  was  equipped  with 
a  spin  chute,  had  a  beefed-up  aft  fuselage  structure  to  handle  the  spin  chute 
loads,  and  had  offsetting  ballast  in  the  nose.  Thus  its  inertia  characteris¬ 
tics  are  quite  different  from  the  average  F-4j,  as  shown  in  the  following: 


F-4e  (Ref.  8) 

f-4j 

stores 

empty  tank 

empty  tank 

wing  stores 

pylons  1,2, 8, 9 

pylons  2,8 

W  (lb 

4o,ooo 

37,000 

c.g.  ii  MAC) 

28.1 

29.3 

Ix  (slag-ft2) 

27,500 

23,850 

ly  (slug-ft^) 

157,000 

127,400 

Ig  (slug-ft^) 

180,600 

146,000 

^Xz  (alug-ft2) 

5,500 

2,  210 

This  difference  is  assumed  to  have  minor  influence  on  the  basic  aerodynamic 
stall/separture  characteristics  as  a  function  of  AOA  but  a  significant 
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influence  on  dynamic  response  parameters  such  as  dutch  roll  frequency  and 
response  to  control  inputs.  With  this  in  mind,  ve  will  proceed  with  a  com¬ 
parison  of  our  P-4j  model  against  the  F-4E  flight  test  results. 

Reference  8  indicates  that  the  clean  aircraft  exhibited  lateral/direc¬ 
tional  stability  breakdown  in  the  form  of  a  slightly  divergent  dutch  roll 
(wing  rook)  in  the  vicinity  of  19-22  deg  AOA.  As  AOA  was  further  Increased 
the  motion  progressed  from  primarily  roll  to  yaw  (nose  slice)  in  the  region 
of  22-25.5  deg  AOA.  The  dutch  roll  mode  evaluated  over  the  same  AOA  region 
from  our  6  OOF  frozen  point  model  is  shown  in  Fig.  29*  This  also  indicates 
that  at  zero  sideslip  the  dutch  mil  slowly  becomes  divergent  at  about  19 
deg  AOA.  The  aero  data  plot,  Pig.  25,  shows  that  0^^  passes  through  zero  at 
about  20  deg  AOA,  while  at  this  same  point  (Pig.  lO)  is  still  relatively 
large.  One  would  expect,  then,  that  the  dutch  roll  motion  would  be  primarily 
rolling  motion  (wing  rock) .  By  25  deg  AOA  is  very  small,  while  0^^  Is 
very  large  negative.  Thus,  one  would  expect  the  dutch  roll  mode  to  be  pri¬ 
marily  one  of  yaw  motion  (nose  slice). 

The  root  locus  of  Pig. '29  also  indicates  that  the  dutch  roll  mode  is 
quite  sensitive  to  sideslip  in  this  same  AOA  range.  Therefore,  it  is  neces¬ 
sary  also  to  look  at  time  trace  comparisons  between  flight  test  and  our  com¬ 
plete  6  DOF  nonlinear  aero  model.  Figure  50  is  an  P-ij’E  time  trace  from  Ref. 

8  at  21  deg  AOA  and  with  all  augmentation  off.  (It  also  has  traillng-edge 
flaps  at  half  deflection  but  they  have  a  relatively  small  influence  on  lat- 
eral/direotional  stability.)  The  traces  show  the  dutch  roll  to  be  a  constant- 
amplitude,  limit-cycle-like  oscillation.  A  similar  set  of  traces,  starting 
at  21  deg  AOA  from  our  P-4j  6  DOF  model  shows  an  almost  identical  oscillation 
(Pig.  5l ) .  A  comparison  of  specific  traces  shows  the  following  peak-to-peak 
excursions. 


Flight  Teat 

Simulation 

a 

21  deg 

21  deg 

P 

1 1  deg 

6.3  deg 

tP 

30  deg 

21  deg 

P 

60  deg/sec 

24  deg/sec 

r 

2.5  deg/sec 

1 .4  deg/sec 

period 

3  sec 

6  sec 

maneuver 

wind-up- turn 

straight/level 
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Figure  29.  F“4j  6  DOF  Linearized  Equation;  Lateral/ 
Longitudinal  Root  Migration  with  a  and  3 


FigiU’e  30.  Controls  Fixed  Wing  Rock 
Time  History  (Ref.  8) 


The  absolute  magnitude  of  motions  depends  upon  the  initial  excitation 
and  is  not  particularly  close.  However,  the  ratios  of  the  motions  are  quite 
close  and  again  tend  to  validate  the  aero  model,  Eie  period  of  the  oscilla¬ 
tion  is  quite  different  because  of  the  higher  energy  (dynamic  pressure)  of 
the  flight  test  maneuver. 

The  nonlinear  response  traces  of  Fig.  3l  show  one  additional  factor. 

It  will  be  noted  that  the  longitudinal  traces  show  a  frequency  which  is  pre¬ 
cisely  double  the  dutch  roll  frequency,  while  Fig.  29  predicts  the  longitu¬ 
dinal  short  period  to  be  higiily  damped  and  about  equal  in  frequency  to  the 
dutch  roll.  Other  F-4E  flight  test  time  traces  from  Hef.  8  (e.g.,  Fig.  32) 
reflect  the  same  "frequency  doubling"  noted  in  our  aero  model.  This  "fre¬ 
quency  doubling"  can  only  be  caused  by  the  pitching  moment  due  to  sideslip, 
Cmp;  which  is  quite  strong  and  negative  for  the  F-4  aircraft.  Note  in  Fig. 
32  that  each  peak  in  the  a  trace  coincides  with  a  zero  crossing  in  p.  In 
this  set  of  traces  the  rudder  activity  indicates  that  the  yaw  SA8  was  on. 

A  comparison  set  of  traces  from  the  F-4j  model  with  pitch  and  yaw  SAS  on, 
starting  at  a  trim  of  Oq  “  23  deg,  »  o  is  shown  in  Fig.  33.  Again,  the 
motions  of  Pigs.  32  and  33  are  remarkably  similar. 

A  final  comparison  between  flight  and  simulation  is  shown  in  Fig.  3^1. 
Both  aircraft  depart  from  a  wind-up  turn  to  the  left. 

On  the  basis  of  the  above  it  was  considered  that  the  aero  model  exhibits 
characteristics  adequately  representative  of  the  P-4  aircraft  for  use  in 
analysis  and  simulation  of  stall/departure. 
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11.11 1  ..  lll|,lilMil.||)jgpAIM|||PPKRV«' 


Figure  34.  Departure  During  Left  Wind-Up  Turn 
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AFHNSn  II 
r-lU  LkU  FAOXAaS 


A.  mnioouraioK 

This  appendix  documents  the  mathematical  model  used  in  the  simulation 
of  the  P-i4A  aircraft.  The  model  represents  the  "clean"  (no  slat  or  flap 
deflection)  aircraft  at  low  Mach  nxunber.  All  aerodynamic  characteristics 
are  for  a  wing  sweep  of  22  deg  (Pig.  35)  which  la  maintained  in  flight  by 
the  Mach-Sweep  Programmer  below  M  =  0.7.  Aerodynamic  coefficients  are 
derived  from  wind  tunnel  data,  with  modifications  based  on  parameter  iden¬ 
tification  using  flight  test  data. 

The  appendix  is  organized  as  follows;  the  equations  for  the  aerodyn¬ 
amic  forces  and  moments  are  presented  in  Section  B;  the  origins  of  the 
individual  aerodynamic  coefficients  is  discussed  in  Section  C;  the  coeffi¬ 
cient  "look-up"  tables  are  described  in  Section  the  coefficient  data 
are  listed  in  Table  12  and  plotted  in  Pigs.  56  through  71 • 

The  equations  of  motion  presented  in  Appendix  I,  Section  B,  are  applic¬ 
able  for  this  model, 

B.  ABHODINAMIC  FORCE  AED  MDMEET 

OOBFFIOIENT  E^ATIONS 

The  aerodynamic  model,  Table  10,  is  largely  based  on  an  F-l4  model  used 
on  the  NASA  Langley  Research  Center  Differential  Maneuvering  Simulator  (DMS), 
Ref.  9,  ^todifications  were  made  to  incorporste  additional  wind  tunnel  data 
and  simplify  the  functional  dependence  on  sideslip  angle  where  appropriate. 
The  terms  ^Cmp^Q(a)  were  added  to  allow  "tuning"  of  the  model 

using  flight-derived  trim  data.  Similar  terras  wore  included  in  the  lateral 
equations  to  allow  inclusion  of  nonzero  sidoforce,  yaw,  and  roll  moments  at 
p  a  0.  Such  momenta,  arising  from  aircraft  asymmetries  and  asymmetric  vortex 
shedding,  have  been  suggested  as  contributors  to  departure  problems. 
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TABLE  10.  F-i4  aerodynamic  FORCE  AND  MOMENT  EQLTATIONS 


A«rodyn»mle  Force  Eguttloni 

°I*  “  CLbASIC^*' P)  +  ACi^^(a)  +  lB*pl/55)  +  ACi,j^(o) 


(b,  > 

;  -10“) 

■  -’°0L8,/“)P. 

+  (6,  +  10)0.  (o) 

»tg 

(Bg  4 

:  -lOO) 

,(«)(|6«pl/55)  + 

^06,(0,) 

Pdb,,  (<*)»• 

% 

(»l  > 

-10“) 

“lOCBj^^(ft)  + 

(6«  +  ’0)CBs,2(®) 

(8.  « 

-lOO) 

^  ‘  P»aA8IC^“'P^*  °yp-o(®)  ♦  Py6,(«)B»  +  Oy5j.(«»P)Br  +  ♦  (b/8Vo)[Cyy{o)r  +  Cyp(a)pl 

(b)  Aerodyn«inlo  Mijtiwnt  Eqiuntlom  (Xp,  In  it  8) 

■  ®"nABIC^*'  P)  ■*■  *  C(5‘c«-XMf)/100](Ci,  00*  a  4.  Op  iln  ft)  +  ACi^p(ft)(  |6|p|/55) 

*  ACi«6,(«)tei  +  <  V2Vo)C„,(ft)q 

-  CMj^^(ft)8,  (8,  >-100) 

At^^(ft)  .  -IOCBj^^(a)  +  (B,  +  tO)0B6^^(e)  (6,  <  -io®) 

■  °"BAB1C^“'P5'*  *  ACn5^(u)6,  iln  (8.2  p)  +  [(Xog  ~  W)/’0'>](o/b)Cy  +  Cnj^lajB* 

+  Cn5^{o,P)6r  +  Cn6^p(ft)6,p  +  (b/2Vo)[Cn,(o)r  +  Cnp(ft)p] 

P)  8»  +  C|^^(a,  P)6p  +  C|j^p(<‘)B|p  +  (b/BVo)[C|^(o)r  ♦  C|^(o 

'OTE;  For  conciseness,  &gtab  has  been  reduced  to  B 
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In  the  DMS  model  the  static  yawing  moment  coefficient  was  treated  as 
a  function  of  three  variables;  a,  p,  and  63  (&b  “  ^stab)*  unusual 

influence  of  stabilator  deflection  on  directional  stability  arises  fi'ora 
an  interaction  between  the  engine  inlet  shed  vortex  and  the  vertical  and 
horizontal  stabilizers,  Ref.  lO.  Examination  of  the  data  indicated  that 
this  effect  could  be  modeled  as  a  yawing  moment  increment  which  was  a  func¬ 
tion  of  a,  63,  and  p  added  to  the  static  yawing  moment  Cn;Qy:^^gjQ(<x,p) .  Rather 
than  treat  this  increment  as  a  three  variable  look-up  table,  it  was  approx¬ 
imated  as 

AGfigg(a)6a  sin  (8.2p  ) 

where  AGng  (a)  is  a  one  variable  look-up  table.  The  sin  (8.2p)  factor  in- 
Eur-es  that  for  63  <  0  (TEIJ)  the  effect  is  destabilizing  for  either  positive 
or  negative  sideslip.  The  longitudinal  stabilator  effectiveness  is  nonlinear 
in  63  in  the  DMS  model.  This  effect  was  retained  in  the  STI  model  but  the 
dependence  on  sideslip  was  eliminated. 

0.  DATA  SOUHOES 

The  aerodynamic  data  package  employed  in  the  F-I^i  simulator  was  put  to¬ 
gether  from  several  sources  and  then  modified  somewhat  to  obtain  an  accept¬ 
able  match  between  the  simulation  and  actual  flight  traces.  The  purpose  of 
this  section  is  to  document  the  origins  of  the  coefficients. 

Although  a  n\imber  of  high  angle  of  attack  wind  tunnel  tests  have  been  con¬ 
ducted  for  the  F-i4,  the  following  tests  were  the  primary  sources  for  the  STI 
model: 


15  March-16  April  1971;  NASA  Ames  Research  Center, 

12  ft  pressure  tunnel 

Atigust  1  971  j  NASA  Langley  Research  Center, 

3o'  X  60'  (full  scale)  tunnel 

The  complete  test  reports  were  not  available  for  these  test  programs. 
Selected  ARC  data  were  available  in  Ref.  10  and  selected  LRC  data  in  Refs.  11, 
12,  and  13.  The  ARC  data  la  generally  considered  to  be  the  "best",  primarily 
because  the  Reynolds  number  is  higher)  however,  its  use  was  limited  in  that 
Ref.  10  contained  only  static  roll  and  yaw  coefficients  as  functions  of  a  and 
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p  and  normal  force  and  pitching  moments  as  a  function  of  a  only.  No  dynamic 
derivative  or  control  effectiveness  data  is  available  from  the  ARC  test. 

These  were  obtained  from  the  Ref.  9  DMS  simulation  which  was  based  on  the 
IiRC  30'  X  6o'  data  with  some  modifications  baaed  on  siimilation  pilots'  opin¬ 
ions.  The  DMS  data  appears  to  be  primarily  "raw"  LRC  wind  tunnel  data  as  it 
shows  characteristic  scatter,  and  coefficients  (such  as  Cfj(a,p))  which  are 
not  symmetrical  in  p.  When  this  data  was  used,  coefficient  values  at  posi¬ 
tive  and  negative  p's  were  averaged,  where  appropriate,  to  produce  tables 
symmetrical  in  p.  In  some  cases  the  averaged  curves  were  further  smoothed 
to  reduce  questionable  scatter. 

The  adequacy  of  the  wind  tunnel  derived  aerodynamic  data  was  validated 
by  comparison  of  6  DOF  analytic  models  with  available  F-'lh  flight  test  data. 

On  the  basis  of  these  comparisons  several  aerodynamic  coefficients  were  fur¬ 
ther  modified  as  explained  in  Section  E. 

Plots  of  the  coefficients  are  shown  in  Figs.  36  to  71  •  The  reference 
point  for  all  moment  coefficients  is  on  the  FRL  at  MAC.  The  spoiler 
produces  no  effect  above  a  =>  10  deg  and,  since  primary  interest  is  in  the 
high  a  I'egime,  all  spoiler  functions  are  zero, 

D,  AZHODYNAMIC  OOEFTIdEN!!!  DATA 

For  use  by  the  simulation  program  all  aerodynamic  coefficients  are  repre¬ 
sented  as  piece-wise  linear  approximations  in  "look-up"  tables.  These  tables 
are  functions  either  of  a  only  or  of  a  and  p  (see  Table  1l).  In  either  case, 
the  a  breakpoints  are  every  5  <ieg  from  a  =  0  to  55  d.es»  The  p  breakpoints 
are  every  5  Aeg  from  p  =  0  to  ±20  deg.  A  listing  of  the  computer  data  file 
containing  the  lookup  tables  is  shown  in  Table  12.  The  order  of  the  individ¬ 
ual  coefficient  tables  ju  the  aerodynamic  data  file  is  as  shown  in  Table  11. 
For  finctions  of  a  only,  a  is  incremented  by  column  and  then  row,  as  indicated 
below: 

0,  5  deg,  10  deg,  15  deg,  20  deg 
25  deg,  30  deg,  35  deg,  4o  deg,  45  deg 
50  deg,  55  deg,  0,  0,  0 


For  functions  of  a  and  p,  P  is  incremented  first,  then  a,  as  follows; 

Row 

1  (0,  -20),  (0,  -15),  (0,  ”10) »  (0,  -5),  (0,  0) 

2  (0,  5),  (O,  .10),  (0,  15),  (0,  so),  (5,  -20) 

22  (55,  10),  (55,  15),  (55,  20),  (o),  (o) 

For  convenience  in  creating  computer  mnemonics,  the  control  variable 
8„4.^v  has  been  abbreviated  as  8^  in  the  data  file  and  plots. 


TABIjE  1 1 .  ORGANIZATION  OF  AERODYNAMIC  DATA  PILE 


PABAMKTER 

ruNCTioM  or 

FWMam 

urns 

an  MaMONtc 
AMO  mict 
SINEN8IOM 

NIAQER  OF 
MTA  PTS 

CABO 

(LfflB) 

CODE 

^LBaSlC 

at  B 

— 

CLUB(ia,9) 

106 

BTIIX 

a 

-- 

ooi<BO(  ia) 

ia 

101 

°l6,1 

a 

ci/0Si(ia) 

ia 

IOU 

a 

l/i*i 

aoaaOa) 

IS 

105 

6CIJP 

a 

— 

DCLarda) 

la 

10£ 

°Qusie 

tti  B 

— 

ecMa(ia,9) 

108 

200 

a 

1/4*1 

assaida) 

ia 

aou 

a 

1/4*1 

ocEaa(ia) 

ia 

205 

^otr 

a 

— 

Dsoarda) 

ia 

ao6 

®yBA8ie 

«!  B 

— 

«BA8(12,9) 

106 

'  >00 

^/B-O 

a 

— 

DcraoCia} 

12 

501 

a 

lA*4 

cvrda) 

18 

50J 

Cvr 

a 

l/pa4 

Cffda) 

18 

504 

a 

1/4*1 

«»(ia) 

la 

JOB 

Oygj. 

a 

1/4*1 

noarda) 

ia 

J06 

OiBMlC 

ASi^aO 


aRBAa(is,9) 


Oip  a 

1/rtd 

0R?<ia) 

12 

4oj 

Oip  «■ 

l/Md 

(}nR(  ia) 

12 

4o4 

Ola,  a.  B 

i/d*a 

aBDAd8,9) 

106 

405 

fliBiP  * 

i/4*a 

CBDeF(  12) 

12 

4o£ 

'la-  ^ 

l/4*f 

CBn(l2>9) 

10a 

4o7 

,  -J 

C»<BA0(l2i9) 

108 

500 

“■BASIC  '  ^ 

OCm.o  a 

— 

DRI«0d2) 

12 

501 

Caq  a 

l/r,d 

0I«(12) 

12 

503 

8*B., 

l/d»« 

CMC81(12) 

12 

504 

o*a,s  * 

l/d*a 

0MDB2d2) 

12 

505 

^®“BP  ^ 

DCMaPt  12) 

12 

506 

OnsAalo  *'  ® 

— 

UKaASd2,9) 

106 

600 

OOnpMJ  » 

— 

ooKBOda) 

12 

601 

ootia,  » 

i/d,i 

OONDSda) 

12 

602 

Cnp  a 

l/rvl 

Ct*P(12) 

12 

605 

On,  “ 

l/rad 

CKR(12) 

12 

604 

One, 

l/d«l 

a<DA(l2) 

12 

605 

Onasp  “ 

1/4*1 

CH&SP(12) 

12 

606 

Oner  1  *■>  ® 

1/4*1 

CM»d2,9) 
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607 

TABLE  12.  AERODyWAMIC  DATA  FILE 


0.<il45000E-01 
0. 145DSOOE+0O 
0.5744194E+00 
0.6024052Ef00 
0. 10017SAEH'01 
0.9a0S749E+00 
0, 1276219E+01 
0,12135a6E+01 
0.U03160E+01 
0.1640450E+01 
0,17(4i4289E+0i 
0,179044(iE+01 
0.1852B92E+01 
0. 18407S2E'101 
0.1770326E+01 
0.174ai22E+01 
0.1666409E401 
0,1(6(68984E+01 
0.1476011E+‘01 
0.1546793E+01 
0,13Q5247E+01 
0.1.394S64E401 
0>)0<ii0<^0«r0ti> 
O.fOttOtiOtfO.r 


'0."i470500£+6(j' 
O.1098S0OE+00 
O.A0240a2E+00 
O.S74<il94E+00 
0.10194B4E+01 
0.921BSA9E+00 
O.1278946E+0i 
0.14394aiE+01 
0,  lS69(i73E+01 
0.1(498947E+01 
0.1740003E+01 
0.1052892E+01 
0.179044*E+01 
0.i870S20E+01 
0.J.733O31E+01 
0.1771812E+01 
0.15<67197E+0J, 
0.  U52013E‘f01 
O.14999O0EfOl 
0.15213192+01 
0.1394564E+01 
0.1385247E+01 


0. >0. >0. >0. »0. y _ 

.0185. .0125) .0135 > .0150. .0165> 
,0145. .0130. .0115. .0040. .0080. 

.0070, .OQAS.O. .0. .0.  _  _ 

To^j .  .  6 09 5 .  .  0 i"l  OTi! O 7 0 1 2 o , 
.0125. .0125. .0135. ,0120. .0110. 
.0080, .0055,0. ,0. .0. 

?rr.‘07TC77UTr5TT; - 

0..0..0..0..0.. 


0. .0. ,0. ,0. .0.  . 
0.1700000E«bl 
0.2295000E-01 
0.47fll333E-01 
0.S320354E-01 
0.1624715E+00 
0. 158"J04F,IE-tO0 
0.34S4817E+00 
0.3118891E+00 
0,S772239E+00 
0.7522102E+00 
0.8174002E+00 
0. 1031058E+01 
0.108B497E+01 
0.1294506E+01 
0, 1243747E+01 
O.1403168t'+Ol 
0. 139B217E+01 
0,16B3551E+01 
0.17407741' +01 
0.1049619t:+01 
0.198:S867F.+01 
0. 19B2491E  +  01 


0.1705000F;:-01 

o.2iaooooe:--oi 

0.5320554E-01 
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Figure  43. 


E.  VALIDATION 


This  section  documents  the  validation  of  the  F-l4  high  angle-of-attack 
aero  model  against  flight  test  data.  Basically,  it  was  necessary  to  deter¬ 
mine  whether  the  simulation  model  would  produce  the  high  a  dynamic  phenomena 
of  the  unslatted  F-14A,  specifically: 

•  Mildly  divergent  dutch  roll  in  the  AOA  range 
12  to  22  deg 

•  Wing  rock  when  AOA  is  maintained  in  the  1 2  to 
25  deg  range  at  low  Mach  Number 

•  Roll  reversal  starting  at  a  =  l8  deg  when  dif¬ 
ferential  tail  is  used  to  roll  the  aircraft 

1 .  Validation  Methods 

Two  basic  methods  were  available  for  comparing  the  simulation  model  with 
the  flight  teat  data.  First,  simulation  time  response  traces  were  compared 
directly  to  the  flight  test  traces.  This  approach  had  limited  success,  since 
the  STI  nnpiloted  time-domain  simulation  program  does  not  have  provision  for 
non-zero  acceleration  initial  conditions  (i.e,,  the  aircraft  is  assumed  to 
start  from  a  steady  trim  condition)  and  does  not  have  the  capability  to  re¬ 
produce  the  complex  control  inputs  of  the  flight  traces. 

The  second  method  involved  extracting  values  of  damping  ratio,  t,,  and  fre¬ 
quency,  (D,  by  locally  approximating  oscillations  in  the  flight  traces  with  a 
second-order  linear-system  response.  In  regions  where  the  control  inputs  were 
fixed  or  could  be  related  to  identifiable  feedback  loops,  flight- derived  ^ 
and  CO  values  could  be  compared  to  open-or  closed-loop  roots  obtained  from  small- 
perturbation  linearization  of  the  simulation  model.  This  second  method  was 
particularly  useful,  since  the  lateral  oscillatory  response  was  of  primary 
interest  and  the  linear  response  parameters  ^  and  cu  can  be  directly  related 
to  the  aircraft's  stability  and  control  derivatives  through  literal  approxi¬ 
mate  factors.  Ref.  1 .  The  approximate  factors  were  rewritten  for  body- center¬ 
line-axis  stability  derivatives.  Table  i3,  rather  than  the  usual  stability- axis 
factors,  so  that  the  response  parameters  could  be  related  directly  to  the  aero- 


dynamic  model  body  axis  moment  coefficients.  Thus  the  significant  stability 
and  control  derivatives  for  a  given  response  parameter  could  be  determined 
and  modified  to  match  the  flight  test  result.  Such  parameter  modifications 
made  by  linear  analysis  were  then  checked  by  extracting  ^  and  aj  values  from 
non-linear  simulation  time  responses. 


TABLE  13.  APPROXIMATE  FACTORS  FOR  KEY  LATERAL - 
DIRECTIONAL  PARAMETERS  AT  HIGH  AOA  USING  BODY 
CENTERLINE  AXIS  DERIVATIVES  =  0,Po  *  O) 
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2.  Flight  Test  Conditions  and  Configurations 

F-i4-  flight  test  data  covering  the  Wgh-a  region  was  obtained  from  the 
manufacturer,  Ref.  15.  From  this  data  six  flight  test  runs  (Table  l4)  were 
found  which  were  suitable  for  comparison  with  the  STI  simulation.  Three  dif¬ 
ferent  aircraft  were  involved  in  these  tests.  One  flight,  No.  199,  was  only 
used  for  comparison  with  the  simulation  trim  routine. 

Reproductions  of  the  flight  test  traces  are  shown  in  Pigs.  7^,  75,  78, 

80,  82,  and  Bl*-.  Aircraft  2  and  3  had  flight-test  nose  booms  with  angle-of- 
attack  and  sideslip- angle  sensors.  The  a  traces  for  Flights  230-5  and  230-6 
(Aircraft  3)  are  "raw"  noseboom  a  (denoted  ajigu  on  the  traces).  Noseboom  a 
is  related  to  true  a  by 

“TRUE  “  “NOSEBOOM 


82 


AIBCBAPT  CCHFIGORATiaSS 


Aircraft  IX  had  no  nose  boom  and  hence  no  sideslip  sensor.  Angle  of 
attack  was  measured  by  the  ARI  nose  probe  and  denoted  on  the  traces. 

The  nose-probe  a  is  related  to  true  a  by 

“true  "  .8122  otfjQsgpRoBE  •'797  d.eg 

The  simulation  model  was  designed  to  represent  the  "clean"  (flaps 
and  slats  retracted)  P-l4  at  low  Mach  number  and  forward  wing  sweep  (22 
deg);  Pig.  35.  The  flight  test  runs  used  for  validation  match  these  con¬ 
ditions  as  closely  as  possible,  specifically; 

•  Plaps  and  slats  are  retracted  except  for 
Bflap  =  —2  deg  on  Plight  No.  236-1 

•  Mach  number  below  M  «=  .35 

•  Wing  sweep,  A  =  19-20  dag 

•  Glove  vane  and  speed  brakes  retracted 

•  No  external  stores 

The  maneuvers  were  performed  at  a  1  g,  either  by  maintaining  speed  and 
AOA  constant  and  exciting  lateral  motions  with  a  doublet,  or  in  a  "i  g  stall" 
in  which  angle  of  attack  la  steadily  increased  as  speed  is  decreased.  Longi¬ 
tudinal  acceleration,  Vj,  varied  from  0  to  -0.4  g,  and  large  negative  flight 
path  angles  were  reached  at  higher  angles  of  attack.  Fairly  large  excursions 
in  mean  bank  angle  occurred  in  some  runs,  but  lateral  accelerations  were  gen¬ 
erally  small.  True  AOA  ranged  from  13  deg  to  32  deg,  which  covers  the  range 
of  interest  for  high  AOA  phenomena  on  this  aircraft. 

At  low  speeds  the  lateral  SAS  gains  are  scheduled  with  AOA  (Pigs.  72  and 
73)  such  that; 

•  yaw  SAS  off  a  10° 

•  ^latstdck  P^^^ed  in  10°  <  a  <  20° 

•  Roll  SAS  and  command  20“  <  a  <  3l° 

augmentation  phased  out 
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Actuator 
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Where  a  as  used,  in  Figs.  72  and  73  Is  the  angle  of  attack  as  measured  by  the 
Alii  nose  probe.  According  to  this  schedule  the  yaw  SAS  should  be  off  for  all 
cases.  However,  Flight  No.  243  appears  to  have  the  yaw  SAS  on  as  will  be  dis¬ 
cussed  later.  Plights  with  roll  SAS  and  ARI  on  were  included. 

3t  Flight  Test  Conperlaons  and  Serivatlve 

Modifications  to  Match  "Wing  Rook" 

Flight  199  (Aircraft  No.  2) 

Flight  199^  Fig.  74,  is  a  i-g  stall  maneuver  in  which  angle  of  attack  is 
increased  steadily  from  10  deg  to  almost  4o  deg  with  a  ramp  (TEU)  horizontal 
stabilizer  input.  This  maneuver  was  used  to  check  the  simulation  trim  values. 

At  a  given  speed,  1-g  trim  values  of  a  and  Bgtab  ^^om  the  simulation  were  com¬ 
pared  to  the'  flight  a  and  Bgtab  traces.  Initial  comparison  using  the  original 
wind  tunnel  derived  aerodynamic  model  gave  good  comparison  with  flight  a  but 
the  flight  &atab  traces  indicated  a  more  traillng-edge-up  deflection.  The 
original  almulation  Cm  (a,p)  was  based  on  the  LRC  wind  tunnel  data  (which  dif¬ 
fered  from  the  ARC  data  by  an  essentially  constant  nose  up  increment) .  A 
ACni^Q(a')  increment  of  -.055  'was  added  to  the  simulation  to  bring  the  pitching 
moment  into  line  with  the  ARC  data.  This  change  resulted  in  a  better  match  of 
the  Batat;  trace,  as  shown  by  the  circled  points  in  Fig.  74. 

Flight  243  (Aircraft  No.  3) 

In  the  time  "slice"  of  interest  (t  «  33-50  sec)  AOA  is  maintained  at  l8  ± 

2  deg.  Fig.  75.  In  this  AOA  region  the  gain  schedule  implies  that  the  yaw  SAS 
and  ARI  are  off  and  that  the  roll  SAS  is  on.  However,  comparison  of  the  yaw 
rate  and  rudder  position  traces  indicates  that  the  rudder  is  correlated  with 
yaw  rate  at  approximately  the  magnitude  and  phase  angle  that  would  result  from 
the  yaw  damper.  The  rudder  pedal  and  lateral  stick  are  both  effectively  zeroed, 
and  thus  the  aircraft  response  can  be  considered  the  free  response  of  the  aug¬ 
mented  airframe. 
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Figure  7^.  Flight  199  Traces 
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The  aircraft  exhibits  a  divergent  "wing  rock”  oscillation.  The  fact 
that  the  envelope  appears  to  grow  linearly  with  time  rather  than  exponentially 
indicates  a  somewhat  nonlinear  response;  however,  approximate  values  of  damp¬ 
ing  ratio  and  frequency  can  be  extracted  as; 

f  =  -.05 

a'  =  1.4  rad/aac 

The  dutch  roll  root  obtained  from  linear  analysis  using  the  original 
(wind  tunnel  derived)  aerodynamic  data  is  stable  and  about  15  percent  lower 
in  frequency  than  the  flight  value.  This  is  shown  (both  SAS  off  and  SAS  on) 
in  Pig.  76.  The  body-axis  approximate  factors,  Table  15,  evaluated  with  the 
oniginal  aerodynamic  data  yield 


tan  6„  Lp 

• 

B 

(-  tanOoT^  Olp)'^® 

■ 

1.21  rad/sec 

^d“d 

B 

-2*j-Yy  -  Lp  +  (n].  -  vf;;  jj 

B 

4“  ( .0*^75  +  .:?758  -  .0970  -  ,3510) 

.0863 

These  approximations 

are 

within  4  percent  of  the  values  obtained  by  numer. 

leal  factorization  of  the  characteristic  equation.  It  can  be  seen  that 
either  an  Increase  in  dihedral  effect,  t)r  a  decrease  in  roll  moment 

of  inertia,  1^^,  would  increase  the  dutch  roll  frequency.  A  reduction  In 
the  magnitude  of  the  roll  damping,  C£p,  or  a  more  positive  value  of  the 
cross  coupling  parameter,  would  reduce  the  dutch  roll  damping.  Such 

changes  were  made  iteratively  using  the  approximate  factors  for  guidance 
until  the  dutch  roll  reasonably  approximated  the  observed  frequency  and 
divergence  rate. 
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Figure  ''(6,  Comparison  of  Oacdllation  Root  from 
Clmulailon  with  Plight  Tost  for  Flight  No.  245 
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These  changes  produced  a  small-perturbation  (linear',  dutch  roll  root  at 

ti\il  =  [-  .071,  1.4l]  for  the  bare  airi'rame,  and  [  cwd]  =  [-.04l, 

SAS  on,  which  can  bn  compared  with  the  flight-derived  root  in  Fig.  76. 

To  check  possible  nonlinear  effects  the  airframe  alone  'Igital  simulation 
was  run,  Pig.  77,  with  the  final  modified  parameters.  The  pproximate  ^  and 
ui  in  response  to  s  pulse  rudder  were  extracted  as  for  the  flight  test  data. 
This  resulted  in  [^,m]  =  [-.045,  1*4]  for  the  bare  airframe  Actually  the 

frequency  started  at  about  1 .45  rad/sec  and  decreased  to  about  1 .35  rad/sec 
as  the  amplitude  of  oscillation  increased.  With  the  SAS  on  the  root  moves  to 
=  [“.024,  1  .40]. 


The  Flight  234-1  maneuver  was  performed  at  an  approximately  constant  AOA, 
Oq  =  15  deg,  Fig.  78.  In  the  time  slice  of  interest,  t  =  15  to  28  seconds, 
the  yaw  SAS  is  off  per  the  a  schedule,  as  can  be  confirmed  from  the  rudder 
position  trace.  The  ARI  is  scheduled  at  60  percent  but  has  no  effect  since 
the  lateral  stick  is  fixed.  The  roll  SAS  which  would  normally  be  on  at  this 
AOA  has  been  switched  off,  as  can  be  seen  from  the  roll  aeries  actuator  trace. 
Thus  from  t  =  19  to  L3  seconds  the  response  is  that  of  the  bare  airframe. 

This  shows  an  oscillating  divergence  (dutch  roll)  with  [^,l^]  =  [-.07,  1.4  rad/ 
sec] . 

Calculation  of  the  dutch  roll  root  using  the  original  (wind  tunnel)  aero¬ 
dynamic  data  again  produced  a  stable  mode.  Estimates  of  parameter  changes 
required  to  match  the  flight  test  divergence  indicated  that  a  level  cf  roll 
damping  comparable  to  the  final  Flight  243  match,  Ca_  =  -.08  1 /rad,  was  re- 
quired.  When  this  C_gp  value  was  used  in  the  nonlinear  simulation  of  the  bare 
airframe.  Fig.  79,  the  approximate  damping  ratio  and  frequency  were  found  to 

be  [^,aj]  =  [-.058,  1.43]. 


The  Flight  256-1  maneuver  was  a  1  g  stall  in  which  angle  of  attack,  was 
steadily  increased  from  a  =  13  to  28  deg  in  the  time  slice  examined,  Fig.  OO. 
The  roll  SAS  has  been  switched  off,  as  seen  from  the  roll  series  actuator 
trace  and  the  lateral  stick  is  fixed,  thus  the  differential  tail  is  zeroed. 

The  yaw  SAS  is  also  off  per  the  a  schedule.  Between  t  =  35  and  ^0  sec  the 
lateral  stick  and  rudder  are  used  to  excite  a  rolling  oscillation.  The 
lateral  stick  is  then  centered  and  held  fixed.  A  rudder  pedal  input  is  made 
at  t  =  1('3  sec.  Also  note  that  the  rudder  apparently  was  used  to  counter  some 
other  distui'bance  since  the  aircraft  is  rolling  to  the  right  and  the  rudder 
input  is  trailing-edge  left.  Thus  the  rudder  input  and  other  disturbance  may 
contaminate  the  response  sufficiently  to  alter  the  apparent  t,  and  os  values. 

As  in  the  Flight  243  trace,  the  envelope  of  the  bank-angle  oscillation  appears 
to  grow  linearly  with  time  rather  than  logarithmically,  but  an  estimate  of  the 
response  parameters  gives  *  [-.085,  1  .30  rad/sec] , 

Matching  of  this  flight  condition  was  done  directly  from  nonlinear  simu¬ 
lation  time  responses.  Fig.  81 .  Matching  frequency  required  a  30  percent  in¬ 
crease  In  Ix  (at  this  heavy  weight  the  aircraft  has  approximately  4,000  lb  of 
fuel  in  the  wings)  and  matching  the  divergence  rate  required  reducing  the  roll 
damping  below  the  level  of  the  wind-tunnel  data  in  the  20  <  a  <  30  dog  region 
(again  similar  to  the  reduction  required  at  13  <  a  <  20  deg).  The  final 
simulation  time  response  shows  a  nonlogarithmic  9  envelope,  similar  to  the 
flight  traces,  which  is  possibly  a  nonlinear  effect  of  the  time-varying  a. 
Oscillation  frequency  derived  from  the  simulation  9  trace  varies  from  1 .4o  rad/ 
sec  at  low  amplitude  to  1.1  rad/sec  at  large  amplitude.  Damping  ratio  varies 
from  “-.1 1  to  -.085. 


Flight  230-^ 


In  the  time  slice  examined,  the  angle-of-attack  was  roughly  constant  at 
a  =  29  deg.  Fig.  82.  Although  there  is  considerable  stick  activity,  the  roll 
SAB  and  command  augmentation  is  almost  phased  out  at  this  AOA  and  thus  the 
differential  stabilizer  is  essentially  zeroed.  The  large- amplitude  rudder 
oscillation  is  due  primarily  to  the  lateral  stick  deflection  acting  through 
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Figure  60.  Flight  236-1 
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Figure  8l .  Bare  Airframe  Simulation  Time 
Responses  for  Flight  236-1 


the  stick-to-rudder  interconnect  (SRI).  From  an  examination  of  the  rudder 
and  roll  tracea  it  cannot  be  determined  immediately  if  the  pilot  was  attempt¬ 
ing  to  excite  the  roll  with  the  rudder  or  was  feeding  p,  cp,  or  both  to  rudder 
to  damp  the  oscillation.  The  latter  possibility  can  be  checked  by  consider¬ 
ing  a  system  survey  (root  locus  and  Bode-siggiej  see  Part  I,  Section  II)  for 
cp Bp  loop  closure  as  shown  in  Pig.  83.  It  can  be  seen  that  at  gains  on 
the  order  of  those  observed  in  the  flight  traces,  =  1  deg/deg,  the  dutch 
roll  is  stable  but  a  divergent  first-order  root  appears.  Since  no  aperiodic 
divergence  is  apparent  in  the  flight  traces,  it  would  seem  likely  that  the 
rudder  input  is  an  open-loop  excitation.  From  the  flight  traces,  the  fre¬ 
quency  is  m  =  1.4-1 .6  rad/sec  and  Ifp/Br  (Jm)|  =  .6  to  1.0  (-4.4  dB  to  0  dB). 
This  is  in  reasonable  agreement  with  the  open-loop  Bode  plot  of  Fig.  83  indi¬ 
cating  that  the  rudder  effectiveness  is  reasonable  at  this  angle  of  attack, 
a  “  29  deg. 

Flight  230-6 

From  t  =  20  to  4o  seconds,  AOA  is  increased  rapidly  with  a  ramp  stabila- 
tor  input.  Fig.  84.  The  rudder  and  lateral  stick  are  essentially  zero  during 
most  of  this  time.  There  is  no  evidence  of  a  divergent  lateral  oscillation, 
which  is  in  general  agreement  with  the  simulation  which  indicates  a  stable 
airframe  in  this  angle-of-attack  region  (a  >  25  deg) .  The  simulation  yields 


a 

25 

^•.097 

1 .22  rad/s eo 

50 

+  .244 

.93  rad/sec 

4.  Roll  Reversal  validation 

Figure  85  shows  the  migration  of  the  u^p  root  with  AOA.  It  can  be  seen 
that  between  a  =  15  S'hd  20  deg  becomes  real  with  one  zero  in  the  right 
half-plane.  This  produces  an  initial  open-loop  roll  opposite  to  that  com¬ 
manded  with  lateral  stick.  It  also  results  in  a  first-order  lateral-direc¬ 
tional  divergence  if  the  cp  -*■  Bjj  loop  is  closed,  F'ig.  86.  This  behavior  can 
be  traced  to  the  change  from  proverse  to  adverse  yaw  at  a  =  16  deg  shown  in 
the  Cng  vs.  a  plot  of  Fig.  87. 


(a).  I/deg 
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